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CALCULATION OF TOTAL CROSS SECTIONS 
FOR CHARGE EXCHANGE IN MOLECULAR COLLISIONS 

Introduction 

Several varied activities were completed during the time 
period of this grant. As the needs and interests of the 
associated NASA resea ch personnel changed, the direction 
of the research under this grant also changed. Areas of in- 
vestigation included nitrogen ion-nitrogen molecule collisions; 
molecular collisions with surfaces; molecular identification 
fxom analysis of cracking patterns of selected gases; computer 

modelling of a quadrupole mass spectrometer; study of space 

. . o 

charge m a quadrupole; transmission of the 127 cylindrical 

electrostatic analyzer; and mass spectrometer data deconvo- 
lution. Each of these areas will be briefly summarized below. 
More detailed reports are contained in the various status 
reports submitted during the grant period. 

A list of the publications submitted under this grant 
is attached. Talks given related to this research and relevant 
conferences attended are listed. Books and equipment purchased 
with funds from this grant are also given on an attached 
page. 

Molecular Charge Exchange in Nitrogen 


A theoretical investigation of low energy collisions 


of nitrogen ions with nitrogen molecules was made using the 
multistate impact parameter method of Flannery, Cosby, and 
Moran.* The method attempts to take into account the vibra- 
tional states of the incident ions and of the products , with 
the target neutral molecules assumed to be initially in the 
ground state. Total cross sections for N2 + on were ob- 
tained. Energies were such that electronic and vibrational 
states were calculated, while rotational states were neglected. 

The charge-transfer process in diatomic nitrogen is re- 
sonant and symmetric so that the waverunction may be expanded 
in terms of the gerade and ungerade eigenfunctions for each 
electronic state of the quasimolecular complex formed. In 
order to include nonresonant channels , the wavefunction is 
expanded in terms of the molecular eigenfunctions of the 
unperturbed Hamiltonian for the isolated molecular systems . 

This expansion is substituted into the time-dependent Schro- 
dinger wave equation for the internal coordinates and the 
charge-transfer total cross section derived. The interaction 
potential is assumed to be spherically symmetric and can be 
obtained from the gerade-ungerade splitting of the vibrational 
motion. A Morse function is used for the gerade potential, 
and with the modified Lennard-Jones potential** s ^ gives the 
interaction energy for the gerade state for this reaction. 

The ungerade potential is given by Sato . 6 

Vibrational overlaps are needed in order to obtain the 

transition amplitudes for the charge-transfer cross section. 

1 7 

Flanner et al ’ found that the excitation defects occurred 
in groups or "bands" separated by approximately one vibrational 
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quantum of energy, Therefore the channels could be considered 
to be degenerate within a given band, greatly simplifying 
the calculations involved. The charge transfer cross sections 
were obtained by numerical integration, 

The low-velocity approximation^ can be used for collisions 
in which the incident velocity is small enough so that the 
excitation defect is approximately zero. Some representative 
results of both the low-velocity approximation and the full 
multistate treatment as well as some additional experimental 
data from the literature are presented in Figs. 1 and 2 for 
comparison with the Langley data. Calculations were made 
for the cases where the incident ion, target neutral, and 
products were all in the ground electronic state. The vi- 
brational states for the incident ion were the ground and 
the first three excited states; for the target teutral, the 
ground state; and for the products, the ground and first 
two excited states. 

Gas-Surface Charge Transfer Collisions 

The suitability of the multistate impact parameter 
method to gas-gas/surface and gas-surface charge transfer 
collisions was studied. 8-10 it was concluded that although 
this method had direct application to the molecular nitrogen 
collisions, extension to collisions where surfaces are present 
would involve time and effort greater than could be justified 
at the present time. 


6 




7 


Trace Gas Identification 


Computer programs were constructed to identify trace 
gases from data from a mass spectrometer designed for the 
Space Shuttle. Mass peak data for 149 chemical compounds 
were obtained from the data bank of the Cyphernetics Cor- 
poration of Ann Arbor, Michigan, via an acoustic coupler on 
the telephone line. Cracking patterns for each compound 
were sorted by computer which then was programmed to print 
out peak height, name, total mass r and allowable parts per 
million for those which had fragment peaks at each atomic 
mass number from 1 to 207. This list was then searched 
and data output for those mass numbers for which there was 
only one compound fragment present, only two, etc. This 
type of information can be used in programming the computer 
for the Shuttle mass spectrometer to identify positively and 
uniquely the compounds which may be present in the sampled 
gases. 


The Quadrupole Mass Spectrometer 

In order to investigate molecular collisions experimen- 
tally, it is often necessary to select ions for mass and 
energy before the collision and to analyze the same properties 
after. This may be accomplished by means of a mass spec- 
trometer and/or energy analyzer. A computer model of the 
quadrupole mass spectrometer was used to investigate its 

properties. Fig. 3 gives the geometry of a quadrupole mass 

, . 11-13 

spectrometer. 


Figure 3. Geometry of the Quadrupole Mass Spectrometer. 
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The differential equations of motion (the Mathieu equa- 
tions) in cylindrical coordinates for an ion passing through 
a quadrupole were solved by a computer code using a fourth- 
order Runge-Kutta method based on Fehlberg's formulas. 14 
Fig. 4 shows the region of stability for a quadrupole. The 

operating quantities A and Q are related to the direct currant 

. 11 

and radio frequency voltages on the pole pieces. The path 
of an ion can be followed through the quadrupole to determine 
its position and velocity at any instant by the computer 
program. The position, speed, and phase of the incident ions 
were varied to simulate an ion beam entering the quadrupole, 
and transmission curves were calculated for selected beams. 

The position and velocity can be determined for those ions 
exiting, or the location where the ions struck the quadrupole 
rods for those not transmitted. 

Fringing fields at the entrance and exit of the quad- 
rupole greatly influence the ion path. In an attempt to 
reduce the effect at the entrance, it was proposed to test 
a quadrupole in which the pole pieces were aligned in a geo- 
metry resembling a funnel, with the entrance larger than the 
exit. Computer programs were developed so that path-dependent 
boundary conditions could be specified corresponding to different 
pole separations at the entrance and exit. Transmissions were 
calculated for various geometries of the quo irupola . For a 
20 pi trajectory, resolution of 100, with ions entering in 
an annular ring parallel to the quadrupole axis, it was found 
that increasing the quadrupole entrance radius caused the 
width of the transmission peak to decrease, as shown in 


Mass Scan Line 




0 0,706 Q 

Stability Diagram 

Figure 4. Stability region for the' quadrupole mans spectrometer. 
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Resolution 100 
20 Pi Trajectory 
0,f} - r ~ 0,6 

Entrance R 0 

^ 1,00000 

^ 1.01053 

" _ 1.01654 


Operating Point 

Figure 5. Transmission of a quadrupole mass spectrometer 




Figs. 5 and 6. Similar curves are given in Fig, 7 for a 
resolution of 500, Transmission peaks for a 20 pi trajectory 
with the ions entering in an annular ring parallel to the 
quadrupole axis fo* non- slanted pole pieces are given in 
Fig. 8, 

In order to investigate the space charge buildup inside 
a quadrupole mass spectrometer, the behavior of an electron 
in a quadrupole was investigated, The differential equations 
of motion for an electron in the quadrupole were solved and 
various paths plotted, Calculations were performed for the 
NASA Gas Analysis and Detection System (GADS). It wa.s found 
that a large fraction of off-axis electrons entering the 
quadrupole under the given conditions would strike the rods. 

The time which any charged particle spends in a quadrupole 
mass spectrometer was obtained as a function of mass. Nu- 
merical values were computed for an electron and for singly- 
charged ions of various masses introduced into the GADS 
quadrupole . 

Charged particles passing from the ion source into the 
quadrupole field region in the GADS must pass through small 
circular holes in plates held at various potentials. The 
electric field and the potential for such a geometry were 
determined for points on an axis passing through the center 
of the circular hole. The electric field wind the potential 
were also obtained for points in the plane of the opening 
inside the circle using Laplace's equation, Gauss's Law, 
and Ohm's Law. The resulting computed numerical values were 
compared with experimental data obtained at NASA- Langley . 
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The 127° Cylindrical Electrostatic Analyzer 


The theory of the transmission of particles through 

1 c n. e 

the 127° cylindrical electrostatic analyzer ’ was studied. 

The geometry of such an analyzer is shown in Fig. 9. Ana- 
lytical expressions for the analyzer and low-current selector 
transmissions as functions of energy and two input angles, 
ir. and out of the plane of incidence, were obtained. A 
uniform distribution of particles across the entrance slit 
was assumed to obtain analytical expressions for the response 
functions for a given angle for the 1?7° cylindrical elec- 
trostatic analyzer for analyzer geometries in which the entrance 
slit is smaller than, equal to, or larger than the exit slit. 

The resolution and the magnification of the analyzer con- 
sidered as an electron lens were also obtained. Results of 
this study may be used to aid the choice of analyzer design 
parameters and to study optimization of the analyzer. 

The effective response function for the 127° cylindrical 
electrostatic analyzer was obtained as a function of incident 
particle energy from the response function by assuming a uni- 
form distribution in incident angle. Separate analytical 
expressions are given in Figs. 10-12 for the effective re- 
sponse functions for analyzer geometries in which the entrance 
slit is smaller than, equal to, or larger than the exit slit. 
Output curves for the 127° analyzer used as an energy selec- 
tor (monochromator) and as an energy analyzer were obtained 
for input distributions which were Gaussian in energy. Shifts 
in peak energy of the selector and analyzer output curves 
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Figure 10. The effective response function T(k) versus relative 

energy k. Entrance slit = .01 r , exit slit = .02 r . 
Curve A for = 2°, B for 3°, and C for 4°. ° 



The effective response function T(k) versus relative 

energy k. Entrance slit = exit slit = .01 r 0 . 

Curve A for cL = 2° , B for 3°, and C for 4°. 
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Figure 11. 
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Figure 12. 


The effective response function TOO versus 
energy k. Entrance slit = ,01 r , exit slit 
Curve A for o(. - 2°, B for 3°, and C for 3. 


relative 


Oil 



were calculated as a function of incident angular breadth 
and slit width. Broadening effects and resolutions were 
calculated for several representative cases. Broadening 
curves for the three representative geometries are given in 
Figs, 13-15. Fig. 16 presents broadening curves for an 
analyser with equal entrance and exit slits . 

Deconvolution of Mass Spectrometer Data 

The accuracy of any physical measurement is limited 
by the instruments performing it. Preliminary investigations 
of the mathematical technique of deconvolution for data en- 
hancement were made. One of the advantages of this kind 
of computer treatment of data is that costly modifications 
to instruments are avoided. The deconvolution techniques 
of iteration and a Fourier-transform-related' method were 
studied. Data from a magnetic mass spectrometer was used 
for initial testing. 

Let h represent the observed distribution of data, 
f the true or ideal distribution, and g the apparatus 
or instrument function. These quantities are related accord- 
ing to the expression 

h = f * g 

where the asterisk represents the convolution of f with g: 

cc 

h(x) = f * g = J f(y) g(x-y) dy 

- CO 

Physically this means that the ideal or true data describing 
the physical process, f, is smeared out or broadened by 
the measuring instrument or response function, g, to pro- 
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Figure 13. Broadening curves for entrance slit = .01 r , 
exit slit = .02 r . -&F is the full width at 

half maximum (FWHH) of the input distribution 
in energy which is Gaussian in form, A T is the 
FWHM of the response function, and AH is the 
FWHH of the output energy distribution. Curve 
A for =2° and B for 4°. Curve C shows the 
broadenifig curve obtained by assuming the response 
function to be a Gaussian. 
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Figure 14. Broadening curves for entrance slit = exit slit 
.01 r . 'AF, AT, and All as defined for Fig. 
13. Curve A for = 2° and B for 4°. Curve C 
shows the broadening curve obtained by assuming 
the response function to be a Gaussian. 
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Broadening curves for entrance slit = .01 r * 
exit slit = . 005 r . Af, A,T m-i A h ° c 
defined for Fig. 1§. Curve A' for =" 2°^^: 

obtain^ V CUrV ^ C Sh0WS the broadening curv 
a elhstt n y . asi " 8 ll,! responaa function to 



Broadening curves for equal entrance and exit 
slits for limiting input angle <K ~ 3°. Ar, 
A T, and AH as defined for Fig.°13. Curve A 
for slits of .005 r Q , B for .01 r D , and C for 
.10 r . Curve D shows the broadening ( curve ob 
tainec by assuming the response function to be 
a Gaussian. 



duca the data which is actually observed by the experimenter, 
h, The problem is to remove the instrument effects and obtain 
the ideal f from the measured h. Mathematically this 
corresponds to solving a Fredholm integral equation of the 
first kind with a difference kernel. 

Iterative Method 

Morrison's iterative noise removal or smoothing technique 17 "^ 9 
treats the data by first smoothing it, then iteratively re- 
storing the non-noisy output and the compatible noise. It 
was designed specifically to prepare data for deconvolution. 

The initial smoothing produces data h n : 

h 1 * h * g 

The nth restoration is giver, by 

h n s h r.-l + (h - VP *8 n > 1 


After the data is smoothed and then restored the desired amount. 


van Cittcrt's iterative deconvolution or unfolding method 
may be applied. The first unfolded f is assumed to be 
the same as h; 


. 3 -; 


£ 0 = h 

while the nth unfolding is given by 

i? n ~ ^n-l + ^ “ ^n-1 ** s) 

Fourier Transform Related Method 

The Fourier transform may be used to deconvolve data 

because the convolution in the function domain 

h = f * g 

becomes a simple multiplication in the Fourier transform 
domain according to the convolution theorem . 3 ' 7 “ 2 ^ Let 
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the Fourier transforms ba represented by capital letters, 
i.e», the Fourier transform of f is F, etc. Then the 
equation in the Fourier domain corresponding to the convo- 
lution in the function domain is 

II = F G 

Rearranging this equation shows that F may be obtained by 

F » II / G 

if G is not zero. The inverse Fourier transform is then 

computed to obtain f . This is called the straightforward 

deconvolution solution. For discrete data, if a standard 

24-25 

fast Fourier transform program is used, the number of 

Fourier transform points computed is an integral power of 
two. Because noise predominates at high frequencies, only 
the low frequency values are used to calculate f. The 
higher ones above a cutoff frequency are set to zero. 

The results obtained from the straightforward decon- 
volution described above may be enhanced by the technique of 
function domain fitting . 40 An artificial function, a, 
is constructed from the function f obtained from straight- 
forward deconvolution and its Fourier transform A obtained. 
The Fourier coefficients of A and F are compared. Since 
noise tends to dominate the signal more at high frequencies, 
the high frequency coefficients of F, which are mainly 
determined by noise in the data, are replaced with the high 
frequency coefficients of A, which are not affected by 
noise. 

Mass Spectrometer Data 

Mass spectrometer data is usually presented in the form 
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of mass peaks, that is, the percent transmission of the ions 
plotted versus mass. The peak height or more generally the 
area under the peak is proportional to the relative abun- 
dance of that mass. If there is any overlap of the peaks, 
peak location and area (or height) will be correspondingly 
in error. Therefore it is desirable to unfold the overlapping 
peaks to obtain completely resolved peaks to pinpoint accu- 
rately the peak location for mass determination and peak area 
(or height) for relative abundance information. 

Mass spectrometer data may contain two or more over- 
lapping mass peaks. Preliminary investigation of the reso- 
lution of these smeared peaks was made, with encouraging 
results. This research will be continued under another 
NASA grant. 
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